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The rhodium doping in the LaCoi-zRh^Oa perovskite series (x = 0.02 — 0.5) has been studied by 
X-ray diffraction, electric transport and magnetization measurements, complemented by electronic 
structure GGA+U calculations in supercell for different concentration regimes. No charge transfer 
between Co 3+ and Rh 3+ is evidenced. The diamagnetic ground state of LaCoOa, based on Co 3+ 
in low-spin (LS) state, is disturbed even by a small doping of Rh. The driving force is the elastic 
energy connected with incorporation of a large Rh 3+ cation into the matrix of small LS Co 3+ 
cations, which is relaxed by formation of large Co 3+ in high-spin (HS) state in the next-nearest 
sites to the inserted Rh atom. With increasing temperature, the population of Co 3+ in HS state 
increases through thermal excitation, and a saturated phase is obtained close to room temperature, 
consisting of a nearest-neighbor correlation of small (LS Co 3+ ) and large (HS Co 3+ and LS Rh 3+ ) 
cations in a kind of double perovskite structure. The stabilizing role of elastic and electronic energy 
contributions is demonstrated in supercell calculations for dilute Rh concentration compared to 
other dopants with various trivalent ionic radius. 
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I. INTRODUCTION 

The cobalt perovskites LnCo03 (Ln = La, Y, 
rare-earths) show various anomalous behaviors that are 
associated with changes of the spin state of octahcdrally 
coordinated Co 3+ ions of 3d 6 configuration. According 
to recent interpretation in the frame of LS-LS/HS-IS 
model, the process can be understood local exci- 
tation of Co 3+ ions from the diamagnetic LS (low spin, 
t%) ground state to closely lying paramagnetic HS (high 
spin, t\ g e 2 g ) states, followed at higher temperature by a 
formation of a metallic-like phase of IS character (inter- 
mediate spin, t\ g a* ) - see e.g. [H-Q ■ In LaCoOs these two 
steps are well separated (T magn = 70 K, Ti_m = 535 K). 
The experiments show that the HS population, after a 
steep increase at T magn , tends quickly to a saturation, 
reaching 40-50% in maximum - see e.g. the combined 
analysis of susceptibility and anomalous expansion data 
in [4|. To account for such behavior, two formal ap- 
proaches have been used. In the first one, Kyomen et 
al. treat the excitations as locally independent events 
obeying Boltzmann statistics and relate the actual course 
of the diamagnctic-paramagnctic transition to lattice ef- 
fects. An activation energy increasing progressively with 
HS population is obtained An alternative approach 
is based on original idea of Goodenough that the sta- 
bilization of HS state is conditioned by presence of LS 
states in the neighborhood. As shown by Km'zck et al. 
[J] using a simple probabilistic model, when nearest HS 
neighbors are forbidden, the activation energy found by 
the fit has an opposite trend - it decreases with tempera- 
ture (or equivalently with HS population) leading finally 
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to a crossover. Such behavior finds a strong support in 
GGA+U calculations Q and also the limiting HS popu- 
lation of about 45% can be naturally explained. 

An uncertainty exists also concerning the formation of 
metallic phase at around 535 K. The gradual course of 
this process, which is in striking contrast with standard 
I-M transition, has been modeled as a thermal excita- 
tion to IS states [f|. However, it is worth mentioning 
that the two-step transitions in LaCoOa does not neces- 
sarily require an existence of three close lying ionic con- 
figurations (LS, HS and IS). A reference can be done to 
phenomenological model of Bari and Sivardiere, who in- 
vestigated the activation from the LS ground state to 
excited HS states including the magnetoelastic coupling 
between the two species [§[ . Depending on the model pa- 
rameters, they obtained a succession of three phases in 
coherence with experimental findings - the mixed LS /HS 
states with gradually increasing HS weight at low tem- 
peratures, the disproportionation of HS-poor and HS-rich 
sites in the intermediate range, and a re-entered homoge- 
neous phase of LS/HS admixed states at high tempera- 
tures. Very recently this scenario of spin-state transitions 
in LaCoOa based on two states only was supported in a 
study of a two-orbital Hubbard model with crystal-field 
splitting using DMFT calculation by Kunes and Kfapek 
p|- They demonstrated that the existence of dispro- 
portionated phase and subsequent homogenization com- 
bined with closing of the charge gap can be understood 
considering the on-site interactions of valence electrons 
and their itinerancy. This suggests that the physics of 
LaCoOs is primarily of electronic (fermionic) origin, and 
clastic interactions between different spin-state species 
have just a stabilizing effect. 

The present study is undertaken with an aim to re- 
solve by experimental and theoretical means the role of 
electronic and clastic interactions in the stabilization of 
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various Co 3+ spin states in the LaCoC>3 type related sys- 
tems, in particular in LaCo03-LaRh03 solutions. The 
studies on Rh-containing compounds were not frequent 
in the past, mainly due to high cost of rhodium, which 
limits potential applications. Nevertheless, since Rh 3+ 
of Ad 6 configuration is isoelectronic with Co 3+ as regards 
the valence electrons, and is exclusively in low-spin state, 
the substitution is helpful in the fundamental research of 
cobalt spin states and experienced an increased interest 
in the last years. Complete LaC01_2.RlLj.O3 solid solution 
was recently studied in Ref. [|[ with focus on transport 
properties at high temperature and in Ref. with focus 
on magnetic properties. Influence of small Rh dopin g on 
the spin state transition in LaCoC>3 was presented in pi| . 

Our work revisits this system, focusing to composi- 
tions with small doping of Rh (x = — 0.08), includ- 
ing for completeness also the x = 0.2 and 0.5 samples. 
The experiments are complemented with GGA+U elec- 
tronic structure calculations using the supercell with reg- 
ular orderings of Co and Rh ions corresponding to doping 
x = 1/16 = 0.0625, 1/2 and 15/16=0.9375. The method 
applied allows us to investigate not only various valence 
and spin states in the model structures but also to op- 
timize local bond lengths and angles that stabilize the 
given configuration. The role of temperature is simulated 
by a change of unit cell volume. In agreement with the 
experimental findings, the calculations support the sig- 
nificant HS populations preserved in the LaCoi_~Rh r 03 
systems down to the lowest temperature, and suggest 
that the stabilization of HS Co 3+ states compared to un- 
doped LaCo03 is to much extent due to elastic interac- 
tions. Namely, the low-energy configurations are based 
on a nearest neighbor correlations of LS Co 3+ states of 
smaller size and HS Co 3+ or LS Rh 3+ states of larger size 
in a kind of the double perovskite AiBB'0§ structure. 



II. EXPERIMENTAL 

Polycrystalline samples LaCoi_a;Rh x 03 with x = 0, 
0.02, 0.04, 0.08, 0.2 and 0.5, were prepared according 
to the following procedure. Hot water solution of metal 
ions prepared by decomposition/dissolution of La203, 
RI1CI3.H2O, and Co(N0 3 ) 2 -6H 2 in 30% HN0 3 was 
mixed with in advance prepared citric acid (CA) - di- 
ethylenglycol (EG) water solution. The molar ratio 
CA:EG was 1:1 and the molar ratio (CA+EG):mctal ions 
was also 1:1. Clear, transparent and voluminous aerogel 
was obtained after evaporation of water at 140°C. The 
gel was pulverised in a mortar, dried, pyrolyzed, and or- 
ganic carbon residue were removed by proper heat treat- 
ment at temperatures 180-450°C in a chamber furnace 
under static air atmosphere. After careful mixing and 
grinding, the final crystallization of LaCoi-^RhrOs per- 
ovskite phase has been done by annealing at 1200°C for 
100 hours in air atmosphere. 

X-ray powder diffraction using a Brukcr D8 diffrac- 
tometer (CuKa radiation, SOL-X energy dispersive de- 



tector, 29 range 20° - 150°) was employed to determine 
the phase compositions and lattice and structural param- 
eters. The diffraction patterns were analyzed with the 
Rietveld method using the FULLPROF program (ll| . 

The magnetic properties were measured by means of a 
SQUID magnetometer MPMS-7 (Quantum Design) over 
the temperature range 4.2 — 300 K and a magnetosuscep- 
tometer DSM 10 (Manics) at elevated temperatures up 
to 800 K. 

Thermal conductivity, thermoelectric power and elec- 
trical resistivity were measured using a four-probe 
method with a parallelepiped sample cut from the sin- 
tered pellet. The electrical current density varied de- 
pending on the sample resistivity between 10 _1 A/ cm 2 
(metallic state) and 10~ 7 A/cm 2 (insulating state). The 
measurements were done on sample cooling and warming 
using a close-cycle cryostat working down to 3 K. 



III. METHOD OF CALCULATION 

The calculations were done with the WIEN2k program 
[l2l] . This program is based on the density functional 
theory (DFT) and uses the full-potential linearized aug- 
mented plane wave (FP LAPW) method with the dual 
basis set. The core states were defined as an electronic 
configuration (Kr, 4d 10 ) for La, (Ar, 3d 10 ) for Rh, (Ne, 
3s 2 ) for Co and as (He) for O atoms. The valence states 
included 3p, 3d and 4s orbitals for Co, and 4s, Ap, Ad and 
5s for Rh. 

All calculations were spin-polarized. For the exchange 
correlation potential the GGA form was adopted (13j |. 
The radii of the atomic spheres were taken 2.3 a.u. for 
La, 2.05 a.u. for Rh, 1.95 a.u. for Co and 1.55 a.u. 
for O. To improve the description of Co 3d and Rh Ad 
electrons we used the GGA+U method, which corre- 
sponds to LDA+U method described in [3, HI] with the 
GGA correlation potential instead of LDA. The values of 
Coulomb and exchange parameters were used the same 
as in the previous studies, U = 2.7 cV and J = eV, 
for both Co and Rh 0, 0] ■ The atom positions were op- 
timized by minimization of the calculated forces on the 
nuclei [161 ] - 

The calculations were performed in order to test two 
specific cases - small doping of Rh in LaCo03 (and sim- 
ilarly small doping of Co in LaRhOs) and x = 0.5. For 
the first case, an ordered arrangement of 1 Rh and 15 Co 
atoms was constructed (x = 1/16), simulating a rather 
isolated Rh atom in Co matrix. The model is charac- 
terized by intersected Rh-Co-Co-Co-Rh-. . . chains run- 
ning along three pscudocubic axes, and a body centered 
supercell 2\/2a p x 2\/2a p x Aa p is used in order to in- 
clude the octahedral tilting (a p is a lattice parameter of 
the cubic perovskite cell). An analogical model struc- 
ture for half doping x = 1/2 is of double perovskite 
type with 1:1 rock-salt arrangement of Co and Rh ions in 
the \[2a p x \[2a p x 2a p cell, characterized by intersected 
Rh-Co-Rh-. . . chains. The inversion symmetry operation 
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FIG. 1: (Color online) (a) Unit cell volume per formula unit 
(error bars are smaller than symbol size) dependence on x 
of LaCoi-^RhajOg (o). The straight dashed line is a linear 
interpolation between x = and 1. (b) Average Co,Rh-0 
bond lengths (□) and Co,Rh-0-Co,Rh angles (o)dependence 
on x. Corresponding values determined in Ref. [g] are shown 
for comparison (x). 
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was only retained in both cells. 



IV. RESULTS 



FIG. 2: (Color online) Resistivity, thermoelectric power 
and thermal conductivity dependence on temperature of 
LaCoi-^Rh^Oa- 



The symmetry of the unit cell is changed from rhombo- 
hedral i?3c for small doping x = — 0.08 to orthorhombic 
Pbnm for x = 0.2—1 in agreement with previous studies 
The average bond distance Co,Rh-0 increases and 
average bond angle Co,Rh-0-Co,Rh decreases as more 
rhodium is inserted in the structure, i.e. the octahedral 
tilting is enhanced with increasing doping, see Fig. \\)p. 
These observations are in agreement with larger ionic ra- 
dius of Rh 3+ (0.665 A in six-fold coordination) compared 
to Co 3+ . Actually, the ionic radius of Co 3+ significantly 
depends on its spin state (0.545 A in LS and 0.610 A in 
HS state), nevertheless Rh 3+ is larger than Co 3+ in any 
spin state. 

The structural evolution of Co,Rh06 octahedra is man- 
ifested in the increase of cell volume with Rh doping, see 
Fig. [T]a,. The dependence shows a positive deviation from 
linearity in agreement with previous studies @ . This 
deviation could also be detected in the dependence of 
bond distances and angles, although it is mostly hidden 
within error bars because the accuracy of these values is 
much lower, since X-ray diffraction is very accurate in 
determination of lattice parameters, but is not as much 
sensitive to oxygen positions. 

Electric resistivity has semiconducting character at low 
temperature for all x, see Fig. (5^. A transition to a more 



conducting state is observed above 400 K for small Rh 
dopings, manifested by a peak in apparent activation en- 
ergy around 500 K, see inset of Fig. [2k,. For x = 0.5, only 
a broad hump is observed in the activation energy. High 
thermoelectric power at low temperature indicates a low 
concentration of carriers, see Fig. [^Jd. The absolute value 
of Seebeck coefficient ranges from 400 to 800 /iV/K. A 
decrease to a metallic-like value ~ 30 /xV/K is observed 
above 400 K for small Rh doping. In contrast, the See- 
beck coefficient for x = 0.5 saturates at a higher value 
of ~ 100 /iV/K. Thermal conductivity is decreasing and 
the phononic peak at low temperature is suppressed with 
increasing Rh content, see Fig. [2J;. 

Magnetic susceptibility measured in DC field of 10 kOe 
is displayed in Fig. [3] The low-temperature transition 
to LS state is shifted down in temperature by approx. 
35 K for x = 0.02 and it is suppressed for x = 0.04. 
It means that inserting Rh into LaCoC>3 destabilizes the 
purely LS ground state of Co. As can be deduced from 
the slope of inverse magnetic susceptibility at the lowest 
temperatures, a nonzero HS population is present start- 
ing from x = 0.04. The poulation of Co 3+ in HS state is 
estimated between 14 — 18% with maximum for x = 0.2. 
The composition x = 0.2 deserves special attention, since 
the sample shows clear signatures of magnetic ordering 
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FIG. 3: (Color online) Molar magnetic susceptibility and 
inverse susceptibility of LaCoi-iRh^Oa. (The spurious in- 
crease of susceptibility at low temperature for x = and 0.02 
is a Curie term due to minor impurity.) The inset shows cal- 
culated populations of HS and LS states. 



TABLE I: The results of the magnetic susceptibility analysis 
in the middle temperature range (150 — 300 K): effective mag- 
netic moment /i e //, Weiss 9, and calculated ratio of the Co 
cations in HS state (phs, S = 2) and LS state (pls, S = 0). 



x(Rh) 


Me// (MB) 


8(K) 


PHS 


PLS 
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below 10 K, including the finite remanent magnetization 
in measurements of hysteresis loops, in agreement with 
the recent paper of Asai et al. [9| . 

With increasing temperature, the slope of inverse sus- 
ceptibility is changed as a result of the excitation of fur- 
ther HS states and interactions among them. Finally, 
practically linear behavior of inverse susceptibility is ob- 
served within the range 150 — 300 K. At a higher tem- 
perature, the samples x = — 0.2 exhibit a magnetic 
anomaly associated with the I-M transition. Its extent is 
gradually diminished and shifted to slightly higher tem- 
perature with Rh doping, and for x = 0.5 the transition 
completely vanishes, in agreement with the behavior of 
the anomaly in transport data in the inset of Fig. [5^,. 

In order to characterize the room-temperature phases 



of LaC01_2.Rha.O3 in a more quantitative manner, let 
us discuss the results of the Curie- Weiss fit of inverse 
susceptibility within 150 — 300 K, which are summarized 
in Table |U It is seen that the effective magnetic moment 
fJ-eff is enhanced for x = 0.02 compared to x = and 
then it takes an opposite trend and is decreasing with x. 
The number of Co 3+ in HS state {phs, 5 = 2) can be 
calculated from the p e f f presuming Rh 3+ and remaining 
Co 3+ ions in the non-magnetic LS state. The calculated 
Phs is reduced from 0.47 for x = 0.02 to 0.41 for 0.08 
while the number of Co 3+ in LS state pls is practically 
constant around 0.51, as though Rh were predominantly 
replacing Co 3+ in HS state for the small doping range 
of x. On the other hand, the extrapolation to x = 
gives pus = 0.49, which is a somewhat higher value than 
Phs = 0-46 actually determined for x = 0. The origin of 
this discontinuity can be rationalized as follows: 

In the frame of the LS/HS thermal excitation model 
[E &[__]: the population of HS states in pure LaCoOs is 
limited to 50%, since probability of nearest HS neighbors 
is strongly suppressed. Because the excitations to HS 
state are dynamic and only a short-range correlated, the 
actual limit should depend on the size of the correlated 
clusters of alternating HS-LS states (or rather on the 
size of disordered boundaries between the clusters with 
Co ions in LS states), and should be somewhat lower 
than 50%. Therefore p e ff corresponding to pns ~ 0.46 
is observed. We suppose, that the clusters of alternating 
HS-LS states, which are stabilized by dilute Rh dopants 
representing immobile elastic defects, are larger than the 
thermally induced dynamic clusters, hence the sum of 
Rh 3+ and HS Co 3+ (x +phs) may approach closer the 
limit of 50%. 

The decrease of p e f / with higher doping of Rh (x > 
0.08) is not so steep. In particular for x = 0.5, the cal- 
culated pus only decreases to 0.35, while pls is reduced 
down to 0.15. We suppose, that the reason is in cumu- 
lation of large Rh ions, which creates a positive lattice 
strain and thus supports HS state on nearby Co sites. 
This effect of Rh dopants for stabilization of HS state at 
neighboring Co ions prevails for higher x — > 1, whereas 
the previously mentioned effect of Rh ions replacing HS 
Co ions in LS/HS ordered regions and supporting HS 
state on the next neighbor Co sites dominates for lower 
x — >• 0. Nevertheless, in both cases the HS states stabi- 
lized by Rh, unlike those obtained by thermal excitation, 
are preserved down to the lowest temperature. 

The dependence of cell volume on x would be linear if 
the room temperature ratio pns '-Pls remained constant 
at the original x = value for all x. Since this ratio 
becomes higher in the doping region x = 0.2 — 0.5 (and 
presumably also for x > 0.5), the observed cell volume 
show a marked positive deviations, see Fig. QJi. 

The Weiss 8 deduced from the susceptibility data over 
the range 150 — 300 K are negative for all x and the 
absolute values are approximately decreasing with Phs-i 
as expected regarding the antifcrromagnctic interaction 
between Co 3+ in HS state. 
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FIG. 4: (Color online) Examples of Co spin states configura- 
tions considered in the GGA+U calculations. 

The previous GGA+U calculations evidenced, that 
two magnetic phases may exist in LaCoOs perovskite, in 
addition to the non- magnetic LS ground state 0, 0, Ell ■ 
The first one is a result of gradual population of HS states 
conditioned by presence of the LS states at the nearest 
Co sites. Nevertheless, some Co(HS) pairs do exist and 
are responsible for prevailing antifcrromagentic interac- 
tions. The HS state in LS matrix is further stabilized 
if Co(HS) site is expanded (breathing- type) at the ex- 
pense of neighboring Co(LS) sites. The mixed LS/HS 
phase tends to a short-range ordered 1:1 arrangement, 
whose long-range formation is presumably prevented by 
the entropy factor. The second phase consists of large IS 
clusters, which may exist within the LS/HS phase and 
finally tend to a formation of uniform itinerant IS phase 
with ferromagnetic coupling. The cell volume expansion 
destabilize the LS ground state and for certain critical 
volumes the energy of LS/HS phase or IS phase become 
lower. 

In the present GGA+U calculation, several spin 
states configurations of Co were tested using the 
Lai6Coi5Rh04s supcrcell (x = 0.0625), namely: 

1. All Co ions in LS state. 

2. All Co ions in IS state. 

3. Co in LS and HS state in 1:1 ratio, with the 6 near- 
est Co around Rh in LS state, see Fig. [4^, displaying 
the nearest neighbors (Co LS) and the next nearest 
neighbors (Co HS). 

4. Co in LS and HS state in 1:1 ratio, with the 6 near- 
est Co around Rh in HS state, see Fig. HJx 

Spin moment of Rh was allowed to vary, nevertheless it 
always converged to LS value. Cell volumes were adopted 
to values expected for various doping by interpolation be- 
tween LaCoOs and LaRh.03. Atom positions determined 
at room temperature were used as the initial values and 
all rcfinable coordinates were optimized. The cell volume 
was also decreased by 1-2% to simulate the cell contrac- 
tion with temperature (T — > K). 

Fig. [5] shows a comparison of the calculated energies 
for the LS/HS configuration and homogenous IS phase 
relative to the LS phase, depending on the perovskite 
cell volumes of LaCoOs and LaCo 15 / 16 Rh 1 / 16 03. 

The cell volume of LaCoOs at 5 K is about 2% smaller 
compared to room temperature fl9j . however only 1% 
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FIG. 5: (Color online) Stabilization energy of LS/HS (full 
symbols) and IS (open symbols) configurations relative to LS 
configuration of LaCoi-^Rh^Os for x = (□) and x = 1/16 
(o) in dependence on cell volume by GGA+U calculation. 



is related to temperature induced contraction and an- 
other 1% is related with the reduction of the Co 3+ size 
at the spin transition [2(j. Since the spin transition 
in LaCoi-ajRh^Os is incomplete, the maximum volume 
contraction is expected between 1 — 2% in this case. 
The relative cell volumes approximately corresponding 
to T = K are indicated in the figure. In both cases, 
the lowest energy state at room temperature (volume 
Vrt) is LS/HS, but the stabilization energy, relative to 
LS state, is higher for Rh-doped compound. With the 
volume contraction simulating the temperature decrease, 
the LaCoOs system undergoes a crossover at about 99.4% 
volume to the LS ground state, while this crossover for 
LaCois/igRhiyigOs is around 98.7% volume, which is al- 
most T — > K in this case. 

These calculations thus evidence that the LS/HS state 
is stabilized upon rhodium doping. The driving force of 
this stabilization is the elastic energy. The ionic radius 
of Co 3+ strongly depends on the spin state. This is con- 
firmed by GGA+U structure optimization, which gives 
Co-0 bond lengths 1.92 A and 1.97 A for LS and HS 
states, respectively. The optimized Rh-0 bond length is 
around 2.02 A. Inserting Rh to Co matrix therefore brings 
about an increase of elastic energy, which is relaxed by 
keeping the nearest Co neighbors in LS state and excit- 
ing the second nearest Co neighbors do HS state, cre- 
ating locally ordered cluster of large cations (Rh 3+ and 
HS Co 3+ ) and small cations (LS Co 3+ ), see Fig. UK- The 
alternative LS+HS state configuration, see Fig.|4}D, is ev- 
idently less favorable regarding the elastic energy, as it 
was confirmed by the calculations. 

Calculation for the x = 0.5 within the La4Co2Rli20i2 
supercell suggests that the LS Co 3+ state is stable for 
T — > 0, energies of HS and LS state become comparable 
at around the room temperature, and HS state is stabi- 
lized at elevated temperatures, i.e. upon further increase 
of the unit cell volume. 

The calculation simulating the rhodium-rich limit x — > 
1, i.e. using the Lai6CoRhis048 supercell (see Fig. 0J:), 
shows that in this case the HS state of Co is more stable 
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FIG. 6: (Color online) Stabilization energy of LS/HS rel- 
ative to LS configuration of LaCoi-zMzOs (M = Co, Rh, 
Al, Ga and In) for x — 1 /16 in dependence on ionic radius 
tm (the average for Co 3+ in LS and HS state is used for 
rco), by GGA+U calculation for room-temperature cell vol- 
ume Vrt(D), V/Vrt = 99% (o) and V/V RT = 98% (A). 



than LS state. This result can be understood considering 
that isolated LS Co 3+ in the lattice of much larger Rh 3+ 
ions would represent an elastic defect of large energy cost. 

The results of GGA+U calculations are in agreement 
with magnetic susceptibility analysis. For small dop- 
ing the Rh atoms stabilize the HS states of cobalt at 
the next-nearest sites due to the elastic energy, therefore 
the low-temperature transition to pure LS state is sup- 
pressed. But at the same time, Rh ions actually occupy 
places which would be otherwise available for HS Co 3+ , 
thus the saturated number of Co in HS state is lowered 
with the doping and number of Co in LS state is retained. 

According to LS-LS/HS-IS model, the high tempera- 
ture spin transition is based on formation of uniform IS 
phase. Thus the destabilization of the IS state by Rh 
doping (Fig. [5]) is in agreement with suppression of this 
transition evidenced by magnetic data. 

In order to investigate the relative role of elastic energy 
associated with larger Rh 3+ size and purely electronic 
effects, additional simulations have been performed for 
x = 1/16 doping of other isovalcnt dopants with differ- 
ent ionic radii - Al 3 + {r Al = 0.535), Ga 3 + (r Ga = 0.62) 
and In 3+ (rj n = 0.80). Cell volumes were adopted to 
expected values for x = 1/16 doping by interpolation be- 
tween LaCoOa and respective LaMOs. Atom positions 
determined at room temperature were used as the initial 
values and all refinable coordinates were optimized. The 
cell volume was also decreased by 1-2% to simulate the 
cell contraction with temperature. The results are pre- 
sented in Fig. [SI It is seen that the stabilization energy of 
the LS/HS configuration with respect to pure LS phase 
is proportional to the ionic radius of the doping cation 
in the Al, Ga and In doped systems, but in the case of 
Rh there is an additional energy gain, obviously due to 
electrons of unfilled shell Ad and their covalency. 

Total and atom projected density of states (DOS) of 
LaCoi-zRhjOs (x = 1/16) for LS spin configuration is 
displayed in Fig. [7J The character of DOS is insulating 



to 
O 

Q 



60 
30 

-30 
-60 
4 
2 

-2 
-4 
4 
2 

-2 
-4 
4 
2 

-2 
-4 





Ai 




total 








1 1 : 


- — " — | 1 1 1 




Co(LS) " 


i , i 








i 1 i 1 


1 1 

Rh 

1 o-C-co^ 


1 — i 1 1 1 — ^ 






1 1 




i.i. 


1 



-4-2 2 

Energy (eV) 

FIG. 7: Density of states for LS Co configuration of 
LaCoi 5 /i 6 Rh 1/16 03. 



60 
30 

-30 
-60 
4 



w 
o 

Q 



tO 

o 

Q 



o 

Q 



CO 

o 

Q 





total 




I 


i 




Co(HS) 

Co(LS) 

mAv 








1 1 1 1 j 


Rh 




— i — i — 1 







I.I. 


1 



-4-2 2 

Energy (eV) 

FIG. 8: (Color online) Density of states for LS/HS Co con- 
figuration of LaCois/ieRhi/igOa. 



K. Knizck et ai, Stabilization of the high-spin state of Co 3+ in LaCoi- x Rh x O^,. 



7 



with a gap about 0.8 eV. The narrow states 0—0.25 eV be- 
low Ep correspond to Rh(t 2g ), whereas Co(i 2 g) and 0(p) 
states are situated lower in energy between 0.25 — 1 eV. 
The states just above Fermi level are mainly of Co char- 
acter. DOS of LS/HS spin configuration is displayed in 
Fig. HI The insulating character of DOS is retained with 
a slightly decreased gap about 0.5 eV. The narrow states 
corresponding to Rh(t2 9 ) are also within the energy range 
— 0.25 eV below Ep, and the states of Co-LS(t 2g ) and 
0(p) could be found lower in energy between 0.25 — 1 eV. 
A sharp peak above Fermi level belongs to Co-HS(t 2g ) 
band. 

It was assumed in the previous discussion, that there is 
no charge transfer between Co 3+ and Rh 3+ . To estimate 
quantitatively charge equilibria in our GGA+U calcula- 
tions we use Atoms In Molecule (AIM) concept of Badcr 
[2l| . In this approach the unit cell is divided into regions 
by surfaces that run through the minima in the charge 
density. The charge on a given site is obtained by inte- 
grating the electronic density within these regions. The 
advantage of this method is that the analysis is based 
solely on the charge density, so it is independent on the 
basis set and atomic spheres used. 

The ionic charges calculated by the AIM method dif- 
fer from the formal valencies due to hybridization be- 
tween cations and oxygen. The ratio of the AIM and 
ideal charge may be regarded as the degree of hybridiza- 
tion. Therefore, a different AIM charge of cations with 
nominally equal valencies can be solely caused by a dif- 
ferent degree of hybridization with oxygen, and cannot 
be directly considered as a charge transfer between the 
cations. Instead, we compare calculated AIM charge of 
Co in undoped and doped structures. The AIM charges 
of cobalt in LaCoOa are within 1.40 — 1.64 depending 
on the its actual spin state, since the hybridization with 
oxygen is predominantly based on e g orbitals and thus 
depends strongly on their occupation, which is different 
for each spin state. Essentially the same AIM charges 
of Co are obtained for doped structure LaCoi-^MzOs 
(a; = 1/16), in spite of the various values calculated for 
the substituting cations M, namely Rh (1.33), Al (1.98), 
Ga (1.94) and In (1.88). It means that the difference 
between AIM charge of nominally Co 3+ and Rh 3+ in 
LaCoi_ K Rh K 03 solely account for the different degree of 
hybridization of Rh and Co, but not present any ground 
for the charge transfer between Co and Rh. 

V. CONCLUSION 

The present experiments supported by GGA+U elec- 
tronic structure calculations provide strong arguments 



that the transition metal ions in the LaCoi-^Rh^Oa 
solid solutions remain in trivalent states. The GGA+U 
calculations suggest that the elastic coupling of the 
nearest and next-nearest cobalt neighbors of the in- 
serted Rh dopant is important in stabilization of the 
spin-state configurations. The Rh 3+ ions are always in 
the non-magnetic LS configuration, while the Co 3+ ions 
may vary, depending on composition and temperature, 
between the non-magnetic LS and paramagnetic HS lo- 
cal configurations. In distinction to pure LaCoOs, in 
which HS Co 3+ states only appear at increased tempera- 
ture through a thermal activation process, the Rh doped 
systems starting from x = 0.04 exhibit certain number 
of stable HS Co 3+ species already in the ground state 
(up to 18% for x = 0.2). The HS population increases 
gradually with increasing temperature and reaches a sat- 
uration above 150 K, similarly to what is observed in 
undoped LaCoOa. The resulting phases, based on the 
LS/HS disproportionatcd cobalt sites, persist at least up 
to room temperature. A gradual transformation to a uni- 
form state with metallic-like conductivity is evidenced at 
elevated temperatures for samples with x < 0.5. 

The main characteristics of the room-temperature 
phase of LaCoOa is the LS/HS population close to the 
1:1 ratio. This can be understood as a result of cooper- 
ative effects, which prefer a regular ordering of smaller 
LS Co 3+ ions and larger HS Co 3+ . Nevertheless, only 
short-range LS /HS correlations are anticipated consider- 
ing the entropy reasons. Similar ordering tendencies also 
exist in the doped systems but the mechanism is differ- 
ent. The primary HS population arises from the elastic 
strain associated with presence of immobile Rh dopants. 
The strain is released by stabilization of LS states on 
the nearest cobalt neighbors and HS states on the next 
nearest ones. These clusters grow as additional HS states 
are thermally activated and tend to a saturation in which 
small LS Co 3+ ions alternate with HS Co 3+ and LS Rh 3+ 
ions of larger size. The stability of such arrangement can 
be related to two factors. The first one is the known prop- 
erty of double perovskite structure A 2 BB'0§ to adapt to 
very different radii of B and B' cations. This is a kind 
of elastic energy optimization, which is effective not only 
for the rhodium but also for other large isovalent dopants 
like In. The second mechanism is of purely electronic ori- 
gin and is associated with incomplete character of the 3d 
and Ad shells of Co and Rh ions. 

Acknowledgments. This work was supported by 
Project No. 202/09/0421 of the Grant Agency of the 
Czech Republic. 



[1] T. Kyomen, Y. Asaka, and M. Itoh, Phys. Rev. B 71, [2] Z. Jirak, J. Hejtmanek, K. Km'zek, and M. Veverka, Phys. 
024418 (2005). Rev. B 78, 014432 (2008). 



K. Knizck et ai, Stabilization of the high-spin state of Co 3+ in LaCoi- x Rh x O^, 



[3] K. Kmzek, Z. Jirak, J. Hejtmanek, P. Novak, and Wei 

Ku, Phys. Rev. B 79, 014430 (2009). [13] 
[4] K. Kmzek, Z. Jirak, J. Hejtmanek, and P. Novak, J. 

Phys.-Condens. Matter 18, 3285 (2006). [14] 
[5] K. Kmzek, Z. Jirak, J. Hejtmanek, P. Henry, and 

G. Andre, J. Appl. Phys. 103,' 07B703 (2008). 
[6] R. A. Bari and J. Sivardiere, Phys. Rev. B 5, 4466 (1972). [15] 
[7] J. Kunes and V. Kfapek, Phys. Rev. Lett. 106, 256401 

(2011). [16] 
[8] J. Li, A. E. Smith, K.-S. Kwong, C. Powell, A. W. Sleight, 

and M. A. Subramanian, J. Solid State Chem 183, 1388 [17] 

(2010). 

[9] S. Asai, N. Furuta, Y. Yasui, and I. Terasaki, J. Phys. [18] 
Soc. Japan 80, 104705 (2011). 
[10] T. Kyomen, Y. Asaka, and M. Itoh, Phys. Rev. B 67, [19] 
144424 (2003). 

[11] J. Rodriguez-Carvajal, P hysica B 192, 55 (1993), [20] 
|http://ill.eu/sites/fullp"rof| (ILL JRC, Version 4.80 - 
Jan2010). 

[12] P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, [21] 
and J. Luitz, WIEN2k, An Augmented Plane Wave + 
Local Orbitals Program for Calculating Crystal Properties 



(Technische Universitat, Wien, 2001). 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 

V. I. Anisimov, I. V. Solovyev, M. A. Korotin, 
M. T. Czyzyk, and G. A. Sawatzky, Phys. Rev. B 48, 
16929 (1993). 

A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. 
Rev. B 52, R5467 (1995). 

F. Iran, J. Kunes, P. Novak, P. Blaha, L. D. Marks, and 
K. Schwarz, Comput. Phys. Commun. 179, 784 (2008). 
K. Kmzek, Z. Jirak, J. Hejtmanek, and P. Novak, J. 
Magn. Magn. Mat. 322, 1221 (2010). 
M. Zhuang, W. Zhang, and N. Ming, Phys. Rev. B 57, 
10705 (1998). 

P. G. Radaelli and S.-W. Cheong, Phys. Rev. B 66, 
094408 (2002). 

K. Kmzek, Z. Jirak, J. Hejtmanek, M. Veverka, 
M. Marysko, G. Maris, and T. T. M. Palstra, Eur. Phys. 
J. B 47, 213 (2005). 

R. F. W. Bader, Atoms in Molecules-A Quantum Theory 
(Oxford University Press, Oxford, 1990). 



